Coupling and coherent electrical control of two dopants in a silicon nanowire 
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Electric control of individual atomd^^ or 
molecules^ in a solid-state system offers a promis- 
ing way to bring quantum mechanical function- 
alities into electronics. This idea has recently 
come into the reach of the established domain of 
silicon technology, leading to the realization of 
single-atom transistors 5 8 and to the first mea- 
surements of electron spin dynamics in single 
donors 9 . Here we show that we can electrically 
couple two donors embedded in a multi-gate sil- 
icon transistor, and induce coherent oscillations 
in their charge states by means of microwave sig- 
nals. We measure single-electron tunneling across 
the two donors, which reveals their energy spec- 
trum. The lowest energy states, corresponding 
to a single electron located on either of the two 
donors, form a two-level system (TLS) well sepa- 
rated from all other electronic levels. Gigahertz 
driving of this TLS results in a quantum inter- 
ference pattern associated with the absorption or 
the stimulated emission of up to ten microwave 
photons. We estimate a charge dephasing time 
of 0.3 nanoseconds, consistent with other types 
of charge quantum bits^H^. Here, however, the 
relatively short coherence time can be counterbal- 
anced by fast operation signals (in principle up to 
1 terahertz) as allowed by the large empty energy 
window separating ground and excited states in 
donor atoms. The demonstrated coherent cou- 
pling of two donors constitutes an essential step 
towards donor-based quantum computing devices 
in silicon. 

Donor atoms in silicon (e.g. P or As) draw increasing 
attention in view of their use as qubits for the devel- 
opment of quant um computing devices based on silicon 
technology^^. Quantum information can be encoded 
either in the charge state of an electron confined in the 
double-well potential of two nearby donors^, or it can 
be encoded in the electron 2 or the nuclear spin state 1 
of a single donor. While charge qubits allow for sim- 
pler and faster schemes for qubit manipulation and read- 
out, spin qubits offer longer quantum coherence^. Re- 
alizing and controlling the tunnel coupling between two 
donor atoms is a basic requirement not only to create a 
double-donor charge qubit, but also to control entangle- 
ment of adjacent single-donor spin qubits^. In the lat- 
ter case, entanglement is provided by a tunnel- mediated 
exchange interaction between the spins of two donor elec- 
trons. Here we demonstrate coherent electrical control of 
double-donor systems in silicon and use Landau- Zener- 
Stuckelberg interferometry±£J to probe the coherence of 




FIG. 1: Device layout, electrical wiring and measure- 
ment scheme. The device consists of a 20-nm-thick, 60- 
nm-wide silicon nanowire etched out of a silicon-on-insulator 
substrate and partially overlapped in its center by two facing 
polysilicon top gates (in gray). The silicon substrate acts as 
an additional back gate. The silicon nanowire is implanted 
with As donors (red spheres) creating degenerately doped 
source and drain leads, while the gate electrodes and the sur- 
rounding insulating spacer layers (green) protect the 40 nm 
long channel region from high-dose implantation. Different 
methods are used to implant a few donor atoms (As or P) in 
the channel region (see Supplementary Material). Transport 
through these few dopants is accomplished by independently 
tuning their levels into the energy window set by the source- 
drain bias (Vsd)- This tuning is accomplished by means of 
the two top gate voltages (Vgl, Vgr) and a back-gate voltage 
(Vbg). A DC current Isd through the dopants is measured. 
For experiments presented in Fig. 3 and Fig. 4 an additional 
microwave signal is applied to one of the top gates. 



charge oscillations between two donors. 

Experiments were carried out on silicon nanowire tran- 
sistors provided with multiple gates as shown in Fig.l. 
Doping impurities were introduced in the silicon chan- 
nel by conventional ion implantation techniques. This 
process was tuned to yield a weakly doped transis- 
tor channel, containing just a few dopants, sandwiched 
between degenerately doped (i.e. metal-like) silicon 
leads^ (see Supplementary Material). We performed 
low-temperature (15 mK) transport measurements in dif- 
ferent samples implanted with either As or P donors, 
as well as in control samples with no intentional chan- 
nel doping. Only in channel-implanted devices could 
tunnel transport through channel donors be achieved. 
The multi-gate device geometry allowed the selection of 
the conductive path corresponding to tunneling through 



2 




pairs of donor atoms connected in series^. This trans- 
port regime was achieved by tuning the transistor into 
its pinch-off state, close to the onset of free-carrier con- 
duction (see Supplementary Material). 

Fig. 2a shows the source-drain current, Isd, of an 
As-implanted device measured as a function of the two 
top gates voltages, Vql and Vqr, for a source-drain bias 
Vsd — 15 mV. Current flow occurs inside two overlap- 
ping triangular regions with the same size and orienta- 
tion, a clear signature^ of transport through two dopants 
in series. To facilitate the analysis of the observed cur- 
rent features we label the donor atom closer to gate Gl 
(resp. Gr) as (resp. Dr) . We denote the corre- 
sponding charge states (Nl,Nr) where Nl (Nr) is the 
electron occupancy of Dl (Dr). In this notation, cur- 
rent flow within the lower triangle of Fig. 2a is accom- 
plished through the (N L ,N R ) charge cycle: (0,0) -» (1,0) 
— >• (0,1) — >• (0,0). Instead, the upper triangle is associ- 
ated with the (Nl,Nr) charge cycle: (1,1) -> (1,0) -» 
(0,1) — » (1,1). Each of the above cycles results in the 
transfer of one electron from source to drain 18 . The off- 
set between the corresponding triangles reflects the en- 
ergy cost of the Coulomb repulsion between the donor 
electrons. We measure an electrostatic energy difference 
of 4.5 meV between the (1,1) and the (0,0) state corre- 
sponding to a dopant separation^ of 30 nm. 

The energy positions of the donor ground-state levels 
are independently controlled by Vql and Vqr- Their 
alignment, schematically shown in Fig. 2b, results in 
a resonant current ridge defining the overlapping bases 
of the two triangles in Fig. 2a (green cross). Moving 
along this ridge corresponds to shifting both levels across 
the bias window while maintaining their reciprocal align- 
ment. Moving away from the ridge results in a finite level 
detuning, eo, between the ground-state levels. 

An additional current ridge can be identified in Fig. 
2a (orange cross) which we ascribe to a resonant tun- 
neling process from the ground-state of to the first 
excited state of Dr (Fig. 2c). From the distance be- 
tween this ridge and the base line we extract an exci- 
tation energy of 7.4±0.4 meV. Due to interface effects^, 
this value is lower than expected for an As dopant in bulk 
silicon 21 , but still an order of magnitude larger than typ- 
ical level spacings in semiconductor quantum dots. In 
addition, current between the ridges is mostly below the 
noise level, meaning that inelastic tunneling between the 
donor ground states is virtually absent. These remark- 
able features indicate that the donor ground states form 
a TLS very well isolated from its environment. 

To further uncover the quantum properties of this TLS 
we shall now investigate the effect of microwave- driven 
charge oscillations between the dopant ground states. 
The hybridization of their ground states, with tunnel cou- 
pling A results in a level anticrossing at eo =0. The en- 
ergy splitting of the lowest double-impurity levels is given 
by \J A 2 + 6q. In the limit of large detuning (|eo| >> A) 
it reduces to eo, and the dopant states are only marginally 
hybridized. The shared electron is strongly localized on 
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FIG. 2: DC transport measurement through two 
As donors connected in series. a Current map, 
Isd(Vgl,Vgr), for a fixed Vsd = 15raV. At 15 mK, the 
observed current features, are confined within two triangular- 
shaped regions, denoting the sequential tunneling of single 
electrons through a pair of As donors in series, labeled as 
and Dr. Vgl and Vgr allow for a full control of the dopant 
energy levels, eo defines the level detuning between the donor 
ground- states, eo = all along the bases of the triangles 
(green cross) where a current peak is observed due to reso- 
nant tunneling via the donor ground states, eo = eVsD at the 
apex of each triangle, providing the scaling factor between 
gate voltage and energy. A second current ridge is observed 
parallel to the base line (orange cross) . This resonance is due 
to the alignement of Dl's ground state with Drs first ex- 
cited state. Its distance from the base line corresponds to 
eo = 7.4 ±0.4 meV yielding a direct measurement of the ex- 
citation energy for donor Dr. We note that no measurable 
current is detected between ground- and excited- state lines, 
indicating a minor contribution of inelastic tunneling between 
the donor ground states, b and c illustrate the qualitative en- 
ergy diagrams corresponding to the green and orange crosses 
in a. 



either Dl or Dr. The corresponding states, \L) and \R), 
define base states of a charge qubit. Coherent oscilla- 
tions between \L) and \R) can be promoted with the aid 
of microwave photons generated by a gate voltage mod- 
ulation in the gigahertz frequency range 22 . These oscil- 
lations result in a measurable stationary current through 
the double donor. 

This effect is shown in Fig. 3b, where Isd is plot- 
ted as a function of microwave amplitude, Vmw, and 
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FIG. 3: Landau-Zener-Stuckelberg interference pat- 
tern under microwave irradiation, a, Time evolution, 
under microwave driving, of the two-level system (TLS) 
formed by the donor ground states. \L) and \R) refer to the 
uncoupled donor ground states and A to their tunnel cou- 
pling amplitude. At eo = the TLS eigenstates are even 
and odd combinations of \L) and \R); for |eo| >> A, they 
localize to \L) and \R). At every passage through eo = 
an "incoming" state \L) (or \R)) splits into a superposition of 
"outcoming" \L) and \R) states. Staying in the same state re- 
quires a Landau-Zener transition across the gap A. Multiple 
passages done within the dephasing time interfere with each 
other, b, Measured Isd(Vmw,£o) interference pattern for a 
microwave frequency Jmw = 15 GHz and a source-drain bias 
Vsd = 5 mV. Vmw, the amplitude of the voltage modulation 
between the two donors, is proportional to the externally ap- 
plied microwave voltage (see below). The static detuning eo 
depends linearly on the position in the (Vgl,Vgr) plane. The 
scaling factor between gate- volt age and energy is adjusted to 
have the horizonal current ridges spaced by the photon en- 
ergy. This factor agrees within 2% with the one detetermined 
from the triangles in Fig. 2a. c, Current peak integrals ex- 
tracted from b. Upward (downward) triangles refer to n > 
(n < 0). The ensemble of 17 data sets is simultaneously fitted 
(black curves) using equation 1 integrated over peaks of order 
n, with three free parameters: A, 77, and the scaling factor for 
Vmw- 



level detuning eo, for a microwave frequency Jmw = 15 
GHz and Vsd = 5 mV. The ground state levels are posi- 



tioned far away from Fermi levels of the leads. The first 
noticeable feature is a set of equally spaced, horizontal 
current ridges positioned at eo = nhfMWi where n = 0, 
±1, ±2,... and h is the Planck constant^. These ridges 
reflect tunneling currents assisted by the emission (for 
eo > 0) or the absorption (for eo < 0) of \n\ photons. We 
have observed well-defined ridges for up to \n\ = 10 at 
Smw = 10GHz (see Fig. 4). 

The second remarkable feature is a strong current mod- 
ulation along each ridge. The resulting triangular-shaped 
pattern of current fringes is a quantum interference effect 
that can be explained as follows. The microwave field 
leads to an oscillatory time dependence of the level de- 
tuning, i.e. io(t) = eo + eVMwcos(27rfMwt), where e is 
the electron charge. For Vmw > e o (condition defining 
the triangular outline in Fig. 3b) the microwave field 
is large enough to drive the TLS through its level anti- 
crossing (at eo = 0) twice per microwave period 1//mw- 
The relative time evolution of the driven energy levels is 
schematically depicted in Fig. 3a. 

After each tunneling event taking an electron from the 
source into Dl, the two-level system is in state \L). At 
the first passage through eo = the system can take two 
possible "paths", ending up either in state \R) or staying 
in the same state \L) (see inset to Fig. 3a). The lat- 
ter possibility corresponds to a Landau-Zener transition 
across the energy gap of the TLS. Quantum mechanics 
allows for these two possible paths to interfere with each 
other at a following passage through the anticrossing. 
This effect, known as Landau-Zener-Stuckelberg inter- 
ference, isanalogous to Mach-Zehnder interferometry for 
photon d 15 l 24 i Each passage through the avoided crossing 
driven by the microwave field is analogous to the passage 
of a photon through an optical beam splitter. During the 
time between two consecutive passages, the out-coming 
"beams", i.e. states \L) and \R), acquire a quantum me- 
chanical phase given by the time integral of the respective 
energies. Their interference will be constructive or de- 
structive depending on whether their phase difference is 
an even or an odd multiple of n, respectively. This phase 
difference is an increasing function of the microwave am- 
plitude, thereby accounting for the observed alternation 
of maxima and minima in the current along each ridge. 

The overall Landau-Zener-Stuckelberg current pattern 
can be modeled by the expression (see Supplementary 
Material): 



r = ^v ^ m 

SD 2 4- /i 2 /4tt 2 T 2 + T 2 (nhf MW - e ) 2 + (1/T R + 1/2T L )A 2 J 2 1 } 



where J n = J n ( e VMW /hfMw) is the n— th order Bessel function of the first kind, T 2 is the dephasing time, 
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fiviw 1 ~ T 2 /3 f MW 1 ~ T 2 f M w 1 ~ 3T 2 
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FIG. 4: Transition from coherent to incoherent driving 
of the double- donor TLS. Current maps for microwave 
driving frequencies from well above to well below the de- 
phasing rate T 2 _1 . a, For Jmw — lOGHz^ 3/T 2 we ob- 
serve sharply defined interference fringes that can be seen as 
the superposition of two interference patterns: one consisting 
of diagonal and anti-diagonal stripes, issued from the inter- 
ference between two consecutive passages through eo = 0, 
and one consisting of horizontal stripes, issued from the con- 
structive interference between consecutive pairs of passages 
through eo = 0. In the latter case, the time delay between 
consecutive pairs equals f M X w an< ^ constructive interference 
is equivalent to having a static detuning equal to an integer 
number of photon quanta, b, For /mw — 3GHz « T^ -1 , co- 
herence is preserved on the time scale of just one microwave 
period. As a result, only the interference stripes due two 
consecutive passages through eo = remain visible, c, For 
fMW = 1GHz w (3T2) -1 , coherence is lost also between con- 
secutive passages through eo = leading to a structureless 
current map. 



Tl (I\r) is the tunnel rate between dopant Dl (Dr) and 
the source (drain) lead. 

At each ridge, Isd(^o) is a Lorentzian function cen- 
tered around eo = nhfMW- Integrating this function for 
each n value yields a set of 17 traces (Fig. 3c) which can 
be simultaneously fitted with only three free parameters: 
A, 77 = T2(1/Tr + 1/2Tl), and the scaling factor be- 
tween the externally applied microwave amplitude and 
Vmw (this factor was already implicitly used to define 
the horizontal scale in Fig. 3b). We obtain A/h = 125 
MHz and 77 = 4.1 ns 2 . 

In order to determine T2, these parameters where fed 
back into expression (1) and used to fit a high-resolution 
measurement of the resonant tunneling current Isd(^o) 
in the absence of microwave excitation. We find T2 = 0.3 
ns, comparable to values reported for other types of 
charge qubits^^. We estimate that dephasing can- 
not be du e to electrical noise through the filtered gate 
leads 25 26 . Instead it could arise from the escape tunnel- 
ing time^l r ^ -1 o r from the noise coming from switching 
offset charges^2Zl. 

Our evaluation of T2 is confirmed by the frequency 
dependence of the Landau-Zener-Stuckelberg pattern^. 
Upon increasing the microwave period f x 



below to well above T 2 , three distinct regimes can be 
identified, as shown by the data sets in Fig. 4. For 
$mw ~ ^2/3 (left panel), multiple passages through 
the level anticrossing occur within the coherence time, 
leading to clearly defined and well-separated interference 
fringes as the ones shown in fig3. For f^ w ~ T2 (middle 
panel), only two consecutive passages are allowed within 
T2 and successive microwave periods are uncorrelated. 
As a result, the interference fringes blur and the nodes 
of the Bessel functions merge into lines departing from 
the zero detuning point. Finally, for f M X w ~ 3T 2 (right 
panel), even consecutive passages become uncorrelated 
leading to an average current with no structure. These 
measurements clearly show the transition from quantum 
coherent to incoherent driving of the double-donor charge 
qubit. 

Because of the large energy separation between the 
donor's ground and excited states, we argue that driv- 
ing signals of much higher- frequency (in principle up to 
fMW — 1 THz, i.e. hundreds of times larger than the de- 
coherence rate) could be used without inducing unwanted 
excitations. Furthermore, recent progress in the position- 
ing of individual donors with nanometer precision^ will 
give access to fine tuning of tunnel coupling between the 
two donors. This high-frequency driving combined with 
a large tunnel coupling should provide the best configu- 
ration for the development of donor-based qubits. 
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